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Abstract 
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1 Introduction 

Among many interesting classical and quantum phenomena in cavities with moving boundaries, the effect 
of photon creation from the vacuum seems to be the most impressive. The first rough (although not always 
quite correct) estimations were made by several authors in Refs. m^lElElElEI- More detailed and consistent 
studies performed for the past decade in the frameworks of different approaches |7ll5ll^ [Tni[m[T^[T^ll4lll5l 
|I3[T7|[T3[Tl|2niiail23E3l21l23ll3l23EH10^ (for other references see review E2) predicted, indeed, 

an exponential growth of the field energy in the case of periodic motion and under the resonance conditions, 
when the walls vibrate at the frequency which is a multiple of the unperturbed field eigenfrequency. 

Now, an experimental verification of theoretical predictions becomes a difficult but realizable task. An 
idea which seems to be the most realistic at the moment is to realize an effective moving mirror by means 
of periodic creation of a conducting layer from a semiconductor film posed on the surface of the cavity wall 
and irradiated by powerful laser impulses Different prototypes of this idea were discussed in [5^05] . 
but not in a periodic regime (for other schemes see, e.g., |36II37| : however, their realizations seem to be more 
difficult). The effective mirror scheme has several advantages over the primary idea of an acoustic excitation 
of harmonical surface vibrations of the wall jj^j . 

The main of them is the removal of severe limitations on the maximal dimensionless amplitude of the 
displacement of the surface, which resulted from tremendous internal stresses inside the wall whose surface 
vibrates at the necessary frequency having an order of several GHz (or higher). In the effective mirror 
scheme, the maximal displacement of the boundary depends of the thickness of the semiconductor film and 
the laser power, and it can be made several orders of magnitude larger than the maximal displacements (not 
exceeding 10~^ cm TS*) which could be achived in the scheme based on the acoustic excitation of the surface 
vibrations. As a consequence, one can diminish (in the same proportion) the number of oscillations of the 
boundary necessary to produce photons in the cavity, thus relaxing the requirements for the Q-factor of the 
cavity and the admissible detuning from the exact resonance ^18 , 19 . Another advantage is a possibility to 
use periodic pulses of arbitrary shape, whose period can be much longer than the period of oscillations of 
the chosen mode of electromagnetic field. 
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However, a practical realization of the scheme depends on its robustness against inevitable experimental 
imperfections, such as, e.g., nonexact periodicity of laser pulses (creating free carriers in the semiconductor 
layer) or fluctuations in their power. Therefore, the aim of our article is to give some estimations of possible 
deviations from ideal results caused by such imperfections. 



2 Basic relations 



The most part of theoretical studies on the nonstationary (dynamical) Casimir effect was performed within 
the framework of the model of a one-dimensional (Fabry-Perot) cavity with equidistant spectrum of unper- 
turbed eigenfrequencies of the field, when all modes are coupled in the resonance case. However, in realistic 
three-dimensional cavities the spectrum is, as a rule, nonequidistant. This results in great simplifications 
in the resonant case, when periodic motion of the boundary excites only one resonant mode, whereas the 
response of all other modes can be neglected in the long-time limit ^3]. Therefore, we confine ourselves 
here to this simplest case (interesting phenomena which could occur in cavities with few accidentally reso- 
nant modes, when the ratio of their unperturbed eigenfrequencies is integral number due to some additional 
symmetry, were considered recently in j^B, 26,29 ). 

In this case, as was shown in the problem of photon generation in the selected mode is reduced to 
the problem of excitation of a quantum oscillator with a time dependent frequency io{t), which is determined 
by the instantaneous geometry of the cavity. For example, in a rectangular cavity with fixed dimensions 
Lx, Ly and a variable dimension Lz(t) we have 
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For fifty years passed after the seminal paper by Husimi (38) . this last problem was studied in detail in 
numerous publications (for reviews see, e.g., (HnEIIl)- It appears that all properties of the quantum oscillator 
are determined, as a matter of fact, by the fundamental set of solutions of the classical equation of motion 



+ u'^{t)x = 0. 



(2) 



Here we actually need only one consequence of the general theory, namely, that the mean number of quanta 
generated from the initial oscillator ground state due to the time dependence of frequency is given by the 
energy reflection coefficient from an effective "potential barrier" given by the function uj'^it). More precisely, 
if uj{t) = LVi for t ~oo and uj{t) = LUf for t oo (this means that the cavity wall is supposed to assume 
some fixed positions before and after the experiment), then one has to calculate the coefficients of the 
asymptotical form of the solutions to Eq. 0, 



x{t ~^oo)= LoJ^^^ [^e"^f' + 7]e-''^f*] 



(3) 



satisfying the initial condition x{t) = lo^ e*'^'* for t 
coefRcients obey the identity 

The mean number of quanta for i — ^ oo can be expressed in terms of ^ and rj as follows |Sll41| 



— oo. Due to the unitarity of evolution, these 

(4) 



M^\r]\'^R/T, (5) 
^ can be interpreted as energy reflection and transmission coefficients 



where R = and T = 1 - i? = |C| 

from the effective potential barrier. 

Eq. shows that no photons can be created during a single pass of the wall from one position to 
another, because for monotonous functions uj{t) the reflection coefficient is limited by the Fresncl formula 



Ri 



(6) 



corresponding to very rapid change of position (for the time much less than the period of field oscillations) 
|42 | . Since real boundaries in laboratory experiments can only move with velocities much less than the speed 
of light, the process is almost adiabatic, and i?i ^ 1. However, it is well known that the reflection coefficient 
can be made very close to unity in the case of periodic variations of parameters due to interference effects. 

For a periodic motion of the wall, the function a;^(i) is also periodic. We assume that it has a form 
like that shown in Fig. The most important assumption is that each "effective potential barrier" is well 
separated from the next one by some interval of time where u = const (i.e., that the wall moves for some 
period of time, returns to its initial position, stays at this position for some time, and then repeats the 
cycle). Thus we exclude monochromatic oscillations of the boundary, which have been already studied in 
detail in 13 . In such a case, each "barrier" can be completely characterized by two complex amplitude 
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Figure 1: A typical time dependence of the effective frequency. 



reflection coefficients and two complex amplitude transmission coefficients, which connect the "plane waves" 
coming from the "left" and from the "right" . Namely, if some "barrier" begins at t = at terminates at 
f = t*, then one can write two independent solutions of Eq. (01 outside the "barrier" as (hereafter we 
consider for simplicity the case of equal initial and final constant values of the frequency w) 

^ (*)-lt(-)e-*, t>U 

X ^-^ut _^^i+)^^^t^ ^^^^ l»J 

(We use the letter t without supescripts to denote the time variable, whereas the same letter supplied with 
supescripts means the amplitude transmission coefficient; we hope this will not lead to a confusion.) The 
coefficients r'^^'> and t^^^ are not independent, because equation (0) is invariant with respect to complex 
conjugation (the frequency uj{t) is real). The following relations hold (see, e.g., I43ll44j : the simplest way to 
obtain these relations is to calculate Wronskians for suitable pairs of independent solutions): 

(9) 

r(-)iW*+r(+)*tW =0, (10) 
|r(-)|2 + |<(-)|2 = |r(+)|2 + |i(+)|2^1. (11) 

Suppose now that we have several well separated "barriers" , the nth "barrier" (possessing coefficients 
ri^"* and ti^^ if it begins at t = 0) being shifted in time by the value i„_i (with = 0) with respect to the 
beginning of the first one. Then one can combine the first n barriers in a single effective barrier, characterized 
by the amplitude reflection and transmission coefficients pi^"* and There are many different methods 

enabhng to express these coefficients in terms of r^^'' and t''j^\ k = 1, . . . , n: see, e.g., ^511^ and references 



therein. We use the method of transfer matrix p iT7III51II51ISnil^l3^IS5| (see also [SJ for historical review 
and [35115^157) for the most recent applications). 

Let us write the solutions to Eq. Q after the fcth barrier (in the region where uj{t) — ui — const) as 



(12) 



with sufhx related to the solution before the first barrier. Evidently, every two sets of the nearest constant 
coefficients, {ak-i,bk-i) and (ak^bk), are related by means of a linear transformation 



bk 



= M, 
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(13) 



Comparing Eqs. p2|) and (|13|l with Q and ©, and taking into account the identities 10- dJ, one can 
express the elements of matrix Mk as follows: 



Mk = 



fk -9l 

-9k fk 



where 



Jk = [tk \ ' 9k = /t^ . 



Consequently, each matrix Mk is unimodular: 

det Mk ^ \fk\^ ^ \gk\^ = I. 

We shall use the notation 

Mn = 



71 ^ n 



(14) 
(15) 
(16) 
(17) 



for the total transfer matrix of n barriers, shifted in time by tfc, fc = 1, . . . , rt — 1, with respect to the initial 
instant t = 0. Applying consecutively the transformations H13|l with account of the phase shifts 6k = Lotk 
(with ^0 =0), one can obtain the following matrix formula: 



Mn = $„-lA/„$,Vi$„_2Af„-l • • ■ $iAf2*lMi, 



where 
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exp(i6'fc) 



Eq. H18() is equivalent to the recursive relations 

Mk = $fc-iMfc$|_^7Ufc_i, 

Fk ^ gkGl_,e^''--^ + fkFk^i, 
Gk ^ gkF^^.e^''"-' + fkGk-i, 
where Fq = 1 and Gq = 0. One can verify that Eqs. (|21|l and H22|l preserve the identities 

\Fk\^ ~ \Gk\' = I. 

The total mean number of quanta (jSJ created after n impulses is nothing but 

The coefficients p„ and t„ themselves obey the nonlinear recurrence relations 

Pn-i + rnSn-i exp(2i6'„_i) 



Pn = 



1 + Pn-lfnSn-l exp(2i6'„_ i ) ' 
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1 + Pn~l1^nSn-l exp(2i6'„_i) 



(18) 
(19) 
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where s„ = r„/T*, pi ^ r[ \ ti ^ t'l 



3 Strictly periodic motion of boundaries 

Formula H18|l becomes especially useful in the case of strictly periodic motion of boundaries, when all one- 
barrier matrices Mk coincide with Mi, and (f>fc = $5^: 

Mn = (27) 

Obviously, the matrix $|Afi is unimodular, as well as matrix Mi. Therefore, we can use the well-known 
formula for the powers of any two-dimensional unimodular matrix S (see, e.g., j48ll54| ) : 

^" = Un-i{z)S - Un-2[z)E, Z = ^TlS, (28) 

where E means the unit matrix and Un{z) is the Tchebyshev polynomial of the second kind. 

In the case discussed, introducing the new parameter ly according to the relations (hereafter 9 = Oi, 
9 = 91, f = /i) 

z = Re {fe-''') = ± cosh i^, (29) 
and using the hyperbolic representation of the Tchebyshev polynomial 

sinh[(n -|- l)i>] 

Un cosh ly) = [\, / ^ , 30 

we arrive at the following explicit expressions for the elements of matrix M.^ (for simplicity, we suppose that 
the sign of variable z in Eq. (|29|) is positive; this sign does not affect the final result \i'6\ for the number of 
created photons): 

G„=,^^e-("-), (31) 

_ sinh(™^ ^e^^n-\) sinh[(n - \)v\ 
smh(t^) sam\y) 

Consequently, the number of created photons after n impulses equals 

M.^\9f-^^. (33) 
smh [y) 

The generation of photons is possible provided parameter v is real, i.e., 

|/||cos(^-e)| > 1, /=|/|e^. (34) 

Since |/| > 1, one has to ajust the phase shift Q to the phase of the inverse transmission coefficient (/?. The 
maximal effect is achieved for 

Q — Ores = (y5 -|- TTm (35) 

(with m an integer). This is equivalent to the following relation between the periodicity of impulses T — /lj, 
the period of the electromagnetic field oscillations in the mode concerned Tq = 2tt/u) and the phase (/?: 

T = ^ (m + ^/^) . (36) 

In particular, for nonmonochromatic oscillations of the boundary, the field mode can be excited not only 
under the condition of the parametric resonance T = To/2, but the period of motion of the boundary 
can be greater than the period of field oscillations. This fact may be important for possible experimental 
arrangements. 

Under the condition H36|l. cosh(i^) = |/| and smh{v) — \g\, so that 

A/'i''"") = s\nh^{nv) (37) 
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where r = r[ ^ is the amphtude reflection coefRcient from each barrier (remember the relations H15|l between 
different coefficients). If |r| <C 1, then one can replace i/ by \r\. Then for n\r\ > 1 we have 

AAr^)«iexp(2n|r|). (38) 

In the case of a harmonically oscillating boundary near parametric resonance, 

u{t) Wo [1 + 2e sin(m)] , n ^ 2(cjo + S) (39) 

{u>o is the unperturbed field eigenfrequency, and Q is the frequency of the wall vibrations), the number of 
quanta created from vacuum is given by the expression |18| (under the conditions \S\ ujq and \e\ <^ 1) 



AA^ ""^ I'-l-j^- (40) 

7 (t^oe) 



For 5 = (exact resonance 15'), (|40|l coincides with H37|l . if one identifies Vf.f = ire (the number of full cycles 
equals n = flt/(2TT) = loqI/t:). 

We see that the concrete form of the law of motion of the boundary L{t) turns out to be not very important 
for the exponential dependence of the number of created quanta on the number of wall's oscillations n for 

^ 1 (although it influences the concrete value of the coefficient ly). What is important, it is the fulfillment 
of the condition (|36|l . 

The generation of quanta in the harmonic case becomes impossible if (5 > uJae. This means that the 
phase deviation from the resonance value accumulated for one cycle, AO — S ■ T ~ 5 ■ tt/loq, should not 
exceed the critical value A^c = its = Vcj- The same result holds, in fact, in the general nonmonochromatic 
case. Indeed, for small values of \g\ <C 1 wc have |/| = + « 1 + |5p/2. Then Eq. (|^ results in 
j/^ « |gp — |A0p, where A9 = 9 — Ores (see Eq. (|35|l for the definition of the resonance phase Ores), so that 



" |gP - 'lAgp ^i'^^' {n^/W^W) ■ (41) 



The critical regular phase shift per cycle equals \AOc\ — \g\ = sinh(j/res) ~ t'res- By the order of magnitude, 
taking into account the Fresnel limit we can estimate the maximal phase deviation per cycle (from the 
resonance value) as \AOc\ < AL/(2L), where AL is the change of the cavity length. 



4 Influence of irregularities of periodic motion 

In real experiments, it can be difficult to maintain the resonance conditions gk = gi, fk = fi, &k = 
kOres exactly. Therefore, it is important to evaluate the influence of systematic and random deviations 
(fluctuations) from these conditions. 

4.1 Amplitude variations 

It seems that fluctuations of absolute values of the effective reflection and transmission coefficients, or 
coefficients \gk\ and \fk\, are less important than fluctuations of phases Ok (i.e., unequal intervals between 
pulses). To show this, let us notice that for identical barriers, the resonance condition H35|) is equivalent 
to the requirement that the second term in the denominators of fractions in Eqs. H25|l and (|26|) is real. In 
such a case, all coefficients p„ have identical phases. Moreover, one can adjust the moments of time tk in 
such a way that this property holds even for barriers of different heights. Under such "generalized resonance 
conditions" , 

Cn+i -Cn- 2ipn + 2{0n+i - On) = mod(27r), (42) 
where r„ — |r„| exp(zCn) and /„ — exp(i(^„), Eq. (|25|l assumes the form 



Writing |r„| = tanh(f„), one can verify that 



\pn\= tanh y^t^fc 



(44) 



(45) 



and the number of quanta in the selected mode monotonously increases with time as 

M - sinh^ . 

For strictly periodic processes (i^„ = v = const) formula H45|) coincides with (|37|l . 
4.2 Phase fluctuations 

Now we suppose that reflection and transmission coefhcients are exactly the same for all barriers, but the 
strict periodicity is broken: 

ek = ke + xk. h = f, 9k = 9, (46) 

where Xk is a stochastic variable. To have a qualitative understanding, what can happen in this situation, 
let us consider first the case, where Xfc is in fact regular function, namely, 

X2fe+1 = +X, X2k = -X, 

which nonetheless destroys a strict periodicity of pulses. This problem still can be solved analytically, because 
each two consecutive barriers can be combined in new effective barriers, forming strictly periodic sequence. 
The total transfer matrix after 2n pulses can be represented in the following form: 
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where 



and / = /e" 
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Using again formula (|28|l . we can express elements of the total transfer matrix 
in terms of the Tchebyshev polynomials of the argument z = Re(/^) + |gp cos(4x). We omit intermediate 
calculations and bring here the final simplified result, which is valid for j^j <C 1, n ^ 1 and |A6'| ^ 1, Under 
these restrictions (which correspond to expected experimental conditions), we can neglect the difference 
between Un{z) and Un-i{z) for z close to 1 and replace cos(2A6') by 1 — 2(A0)^. Moreover, we can neglect 
the diference between the values Af2n and Af2n-i- Finally, we arrive at the following simple formula: 



cos2(2x)sinh2 (|5|nVcos2(2x) - iA0/\g\ 
cos2(2x) - {Ae/\g\r 



(47) 



We see that the presence of phase "jumps" diminishes the rate of increase of the number of photons, 
although the growth remains exponential for sufficiently large number of pulses. The admissible level of 
jumps is determined by the inequality |cos(2x)| > |A0|/|(7|, which shows that "random" deviations from 
the resonance are much less dangerous than systematic ones. Indeed, the maximal admissible value of the 
systematic phase shift per period equals A9max = \g\, in complete acordance with evaluations made at the 
end of section|31 On the other hand, for zero systematic shift from the exact resonance {A6 — 0), exponential 
growth of the number of photons can be observed for any amplitude of "phase jumps" Xi although the growth 
becomes slower with increase of x in the interval (0, 7r/4): 



A/; «sinh2(|5|n|cos(2x)|). 



(48) 



We have performed numerical experiments on the basis of recursive relations (|21|l and (|22|l . choosing 
parameters 9 and / = yl + g^ in Eq. 14t)l) real (so that the resonance condition is Ores — T^rn; obviously, 
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Figure 2: Mean number of photons after n — 500 pulses (with |g| — 0.01) versus the maximal amplitude of 
random phase fluctuations x for fixed values of systematic phase deviations 9. 
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Figure 3: Mean number of photons after n = 500 pulses (with \g\ 
deviations 9, for fixed values of the amplitude of phase fluctuations x- 



0.01) versus the systematic phase 



this choice does not affect the qualitative picture), namely, we chose g = 0.01 (in order to reduce the time 
of calculations and make the results more vizualizable) . The coefficients Xk were chosen with the aid of a 
program generating random numbers in the interval (— XiX): for different maximal values x systematic 
displacements 9 = A9 from 9res = (which is equivalent to 27rrn, of course). The results are given in the 
figures EE 

Fig. 121 shows the mean number of photons in the cavity after a fixed number of n = 500 pulses (which 
generate approximately sinh^(5) ~ 5500 photons for g = 0.01 in the case of strict resonance) versus the 
maximal amplitude of random phase ffuctuations x for fixed values of systematic phase deviations 9. The 
averaging was performed in all cases over 700 realizations consisting of 500 random numbers Xk- The average 
values (xfc) = X]fc=i Xk varied from 10~^ to 10~^ for each realization (so that no systematic shift was 
added due to random choice of Xk)- Fig- El gives the mean number of photons versus 9 for different fixed 
values of x- Fig- Ql shows the relative (and in the insertion - absolute) mean-squared deviations from the 
average numbers, as functions of x for ^ = 0. In Fig.|Sl the upper curve shows the dependence of the critical 
amplitude of random phase ffuctuations (defined by the condition that the average number of created photons 
after n — 500 pulses is between 9 and 10) on the regular phase shift 9; the lower curve gives the maximal 



amplitude of phase "jumps", when the photon generation stops according to Eq. l|T7|l . It is clearly seen, 
that even very high level of random phase fluctuations does not stop the generation of photons, contrary 
to systematic phase shifts. Moreover, the critical amplitude of truely random fluctuations remains nonzero 
even for systematic phase shifts exceeding the value \g\. 




Figure 4: The relative mean-squared deviation of the number of created photons versus the maximal 
amplitude of random phase fluctuations x for fixed values of systematic phase deviations 9 (for ri = 500 
pulses). The absolute mean-squared deviations are shown in the insertion. 
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Figure 5: The dependence of the critical amplitude of random phase fluctuations on the systematic phase 
shift 9 (the upper curve), compared with the maximal amplitude of "phase jumps", when the photon gen- 
eration stops according to Eq. (|47|) (the lower curve). 



4.3 Analogy with resistance of disordered chains 

It is worth noting that the problem of transmission through a sequence of one-dimensional barriers with 
randomly varying parameters (in particular, interbarrier distances) has been extensively studied in the 
context of solid state physics [SH]. There, one of objectives is the resistance TZ of some chain of scatteres, 
expressed (in dimensionless units) by exactly the same formula Q which gives the mean photon number in 



our case. In particular, Landauer [SHI ^sed a consequence of Eq. 



Pn-iP + k«P + 2|r„p„_i| cos(^'„) 



l-|p„P (l-|p„_iP)(l-|r„|2) 

*n = C« + Z^n^i + 2arg(s„_i) - arg(p„_i), 

and, supposing that the phase shift 'i>„ is a random variable, uniformly distributed over a large interval, 
performed averaging over directly in Eq. H49|) . having replaced the term with cos(^'„) by zero. One can 
easily verify that the change of variables, |p„p = tanh(y„), jr^p = tanh(z„), transformes the set of recurrence 
relations H49() with cos(^'„) replaced by zero, to the equations ?/„ = y„-i + z„, whose obvious solution is 
Un = X]fe=i -^k- Thus we arrive at the following formula for the resistance of a chain with completely random 
distances between scatterers, uniformly distributed over sufficiently large intervals: 



fc=l ' ' 

Here |rfcp is the energy reflection coefficient from the fcth obstacle (barrier) and n is the number of obstacles. 
If |r„p — |rp = const, then (|^ goes to the Landauer formula fSSI 

It is worth comparing Eq. H51() with Eq. (|37|l , paying attention to different powers of the reflection coefficient 
\r\. For n|rp > 1 formula (|51|) yields 

A/'i™"'^«^exp(2n|rn. (52) 

Consequently, totally random phase fluctuations without systematic phase shift must also lead to an expo- 
nential growth of the photon number (chain resistance), although with much more slow rate than in the 
resonance case, if \r\ ^ 1. For example, for \r\ — 0.01 one needs about 500 pulses to create 5000 photons 
under the resonance conditions, according to Eq. whereas Eq. suggests that the same number 

of photons can appear after about 5 • 10^ pulses in the case of completely random fluctuations (with zero 
systematic phase shift). 

However, removing the term proportional to cos(^„) from the right-hand side of Eq. 14911 is not well 
justifled operation, because averaging over fluctuations of phase must be performed in the final expression 
for Afn, and not at some intermediate steps (as soon as {f{x)) ^ f{{x))). Moreover, if the random variable 

is uniformly distributed over some interval (— x, x), then the average value (cos(4'„)) = sin(x)/x is equal 
(or close) to zero either if the ratio x/"" is an integer or if x 3> 1. In Fig. Elwe show the results of numerical 
calculations of the dependence of the natural logarithm ln(A/'„) of the average number of photons created 
after n pulses with completely random phase Xfc (defined by Eq. (|46|l '). uniformly distributed over some large 
interval (— x, x)j for x > 1 and 9 — (i.e., in the absence of systematic shifts). For each choice of x, averaging 
was performed over 20 realizations. Qualitatively, we see an agreement with Landauer's predictions: the 
growth of the number of photons (or equivalent resistance) is close to exponential, although the necessary 
number of pulses is much greater than in the case of x < 1 illustrated in the Figs. 12151 However, the results 
of numerical experiments are close to the asymptotical Landauer's dependence H52|) . In Afn = 2n|rp — In 2 
(given by the dashed straight line labeled as "Landauer"), only for integral ratios m — x/t^ (take into account 
the logarithmic scale), whereas the inclinations of curves corresponding to half-integral values of the ratio 
x/tt = m -I- 1/2 (when |(cos(xn))| = [TT{m + 1/2)]^^) are much bigger for small and moderate values of m, 
approximating Landauer's limit only when m > 100. 

The results of this section are in qualitative agreement with other studies, e.g., |31| (where quasiperiodic 
motions of the boundaries were considered), j55| (parametric resonance in periodic paraxial optical systems), 
or |S31|5ni|Hnj (devoted to quantum oscillators with fluctuating frequencies). A possibility of generating 
coherent phonons in solids with time-dependent lattice strains was studied in jfilj . Mathematically, this 
problem is reduced to quantum oscillator with time-dependent (fluctuating) frequency. Another mechanism 
of producing coherent phonons by femtosecond laser pulses, which is mathematically equivalent to the 
problem of quantum oscillator with time-dependent force, was considered recently in 62 . 
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Figure 6: The natural logarithm ln(A/'„) of the average number of photons created after n pulses with 
completely random phase Xk, uniformly distributed in the interval (— x,x), for ^ 0, g = 0.01 and different 
values of x- The straight line with label "Landauer" corresponds to the Landauer's dependence (|52|l with 
\r\ = g- 



5 Conclusion 

We have considered a simple model of photon generation from vacuum due to nonstationary (dynamical) 
Casimir effect in a cavity with noncquidistant spectrum of eigenfrequencies, which is reduced to the problem 
of quantum harmonic oscillator with time-dependent frequency. We have shown that periodic displacements 
of the (effective) boundary must result in exponential growth of the mean number of quanta (photons) 
in the selected mode under certain resonance conditions. The oscillations of the boundary need not to be 
monochromatic, but they must be close to periodic. However, although the period of motion of the boundary, 
Tb, should be adjusted to the period of field oscillations, T/, and to some phase, depending on the concrete 
form of pulse, Tb may be much greater than Tj. This result seems to be very important from the point 
of view of facilitating performing experiments. The concrete form of the trajectory of the wall enters the 
final result only through effective reflection and transmission coefficients from some "barriers" in the time 
dependence of the eigcnfrequency. 

We have studied the influence of deviations from the strict resonance, caused by systematic and random 
perturbations of the optimal resonance trajectory of the boundary. It appears that amplitude perturbations 
are less important than the phase ones. In turn, among the phase perturbations, the most dangerous for 
preserving the regime of photon generation are systematic deviations, which should not exceed rather low 
level of the order of the frequency modulation depth. On the contrary, random phase fluctuations can 
have much greater amplitude (even of the order of unity for zero systematic deviations) without qualitative 
changes in the behaviour of the system. This circumstance also seems to be important for planning future 
experiments, showing that some requirements might be not so rigid as one could suppose, thus diminishing 
the number of technical problems to be solved. 

The results obtained can be easily generalized to the case when the initial state of the field mode was 



not vacuum, but a thermal state. Actually, one should simply multiply the right-hand side of Eq. JSJl by the 
factor 1 + 2{nth), where (nth) is the mean number of thermal photons in the initial state |51 l23[l24ll27) . 

However, the presented study should be considered only as a qualitative model, because many important 
things have not been taken into account. For example, we considered the excitation of a single selected mode, 
neglecting its possible interactions with other modes. Such an approach seems to be justified for cavities with 
nonequidistant spectra, when the boundary performs harmonic oscillations at the resonance frequency 
But the spectrum of anharmonic oscillations contains many frequencies. On the one hand, this fact, perhaps, 
explains why photons can be generated even if the period of displacements of the boundary is greater than 
the period of the field oscillations. On the other hand, in the presence of many frequencies some other modes 
can be excited, too. This question needs a detailed investigation. Also, the effects of polarization of true 
electromagnetic field (not its scalar model) could be important, as was shown in other examples |29II63| . 
And, of course, the problem of consistent account of losses in cavities with moving nonideal mirrors must be 
solved (perhaps, following the lines generalizing the approaches of Refs. |18ll28ll5^ ). 
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